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ON INVARIANCE OF DOMAINS WITH 
SMOOTH BOUNDARIES WITH RESPECT TO 
STOCHASTIC DIFFERENTIAL EQUATIONS 

Vitalii A. Gasanenko 



Abstract. We prove constructible sufficient conditions of lack of exit by 
solutions of stochastic differential Ito's equations from domains with smooth 
boundaries. 



Consider a stochastic differential equation for process G R n . 

n 



fe=i 



here 

a(t,x) := (aj(t, x), 1 < z < n), bk(t,x) := (bki(x), 1 <i < n), x G 

It is assumed that there is such constant L that for functions bij (t, x), di(t, x) 



the following conditions take place 



\a{s,x) - a(s,y)\ + ^ \b k (s,x) - b k (s,y)\ < L\x - y\ 
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n 



a 



(s,x)\ 2 + Y J \h(s,x)\ 2 <L 2 (l + |x| 2 ). 



(2) 



k 



for all x,y G it"\ 

Here | • |, be norm (length) of vector. 

It follows from [ 7, p. 480] that this is sufficient conditions for existence of 
unique solution of (1). 

Let there be given a measurable set K G R n . A set K is said to be invariant 
set of equation (1) , if under condition P(£o £ K) = 1 the following equalities 



The property (3) of trajectories of solution of (1) sometimes is called viability. 
The necessary and sufficient conditions of viability was proved for the first time 
in [1]. Such conditions was proved for more general constructions of equations 
in [2]. The methods of investigations of these articles are different but the 
set K is the same: it is convex and closure. 

This problem was reduced to viability of ordinary differential equations with 
help approximation theorems Ikeda- Nakao- Yamato for homogeneous stochas- 
tic differential equations in [4]. 

The conditions of viability were formulated in terms of asymptotic behavior 
of distance to considered closed set. The analogy conditions of viability were 
proved for inhomogeneous stochastic differential equations and relative closed 
sets in [5]. We observe that test of conditions in terms of distance to sets 
requires the additional investigations. They are checked for convex sets 

effectively. For example, it was done in [5]. 

It was proved necessary and sufficient conditions or only sufficient 
conditions of viadility (3) in [3,6] by probabilistic methods for the specific 
domains K. 

Our purpose is to obtain verifiable sufficient conditions of viability (3) for 
domains with smooth boundaries. Our method of investigation be different 
from other. It is based on the use of Ostrogradskii- Gauss theorem. 

Consider a closed set K in R n with boundary S{ or dK). 



hold 



P(£(t) G K) = 1 for all t > 0. 



(3) 
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Let U (x, r) denote open ball with center in point x and with radius r. The 
union of balls with centers in K is called e - neighborhood K e of the set K : 
K e = IJ U(x,e). We will denote by S e the boundary of K e . 

x<EK 

We introduce the following function 



u e (x) = 



e 2 

Ce g c 2 -m 2 5 if |x| < e, 
0, if |x| > e. 



The constant is choosed c e such that the following equality holds 

J uj e {x)dx = 1. 

Thus 



c e e n j e iH€l*df = 1. 
I€l<i 

If %(•) be characteristic function of set K 2e then for any e > the function 

rje(x) = J x(z)uj e (x - z)dz. 

satisfies the following relations [8, p. 89]: 

< r] e (x) < 1, 7] e (x) = 1, x E K e , 
le (x) = 0, x i K 3e , le (x) E C°°(R n ), \4 a) (*)\ < L a e~W. (4) 



The next statement follows from the axiom of continuity: 

Statement. If C, be random vector in space R n , then the following 
representation takes place 



Er) e (() = P (C e K) + l e , where l e > 0, l e -> 0, when e -> oo. (5) 
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Lemma 1. If P(£(0) G K) = 1 and for some number eo > and any num- 
bers < e < eo the following inequality takes place 

EvMt)) > ^ett(O)), « > 0. 
then the following equality is true 

Proof. Let the condition of Lemma be fulfilled but statement of Lemma don't 
fulfill. If statement of Lemma don't fulfill then there exists such t* that 

P(£{U) G K) < 1 (6). 

Futher according to the statement (5) and the condition of Lemma 1 we 
have the the following inequality in point t* 

P(Z(U) G K) + l 3 , e > P(£(0) G K) + h,e- 

Letting e — > , we arrive at 

P(au) eK)> P(£(0) g K) = 1. 
The latter one contradicts to (6). This contradiction proves the Lemma 1. 

We make the following assumption: the boundary of dK e belongs to class C l , 
when the following condition of smoothness of boundary of K e holds for e < eo 
under small eo > 0. The intersection of boundary of set K e with ball U(x, e), 

x G K 2e ■ 

A e (x) := U{x J e)DdK e 

is surface whose equation in local coordinates (yi, . . . , y n -i) with origin of 
coordinates in point xq G A e (x) has form y n = <p(yi, . . . , y n -i)- 

The function ip belongs to class C l in region D e , which is projection of A e (x) 
on the plane y n — 0. 

Let us denote by v{z) = (vi(z), i = 1,2, ...,n) the unit vector of external 
normal to boundary S in point z G 5. 

It is known, that if surface is given by relation Q(yi, . . . , y n ) = , here Q(-) 
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be smooth function, then the unit vector of normal n has the following form 



n 



Thus, if the if is differentiable ,then the next reprezentation for v{z) takes 
place locally 



\ 



'i+ E ^ .1+ E % 



-. , /c = 1, n — 1 



i<n— 1 



i<n — 1 



/ 



Suppose now that the bondary i^ e , eo > e > under some eo > belongs to 
class C 2 . 

Theorem 1. If the following conditions are fulfilled 

1. The functions a(t,x),bk(t,x),l < k < n, in addition to properties (2), 
under fixed t belong according to classes C 2 (R n ), C 3 (R n ). 

2. sup sup I ^ bji(s, z)vi(z) I = o(e), e — > 0. 1 < j < n, s > 0; 

3. limsup sup ^2ai(s,z)ui(z) - § E ^ fc ^ s,z) Vi(z)b k j(s, z) < 0. 

s>0 zGS e y i i,j,k 3 J 

then (3) takes place. 

Proof. Applying the Ito's formula, we get the following equality 

t 



Here 



n d 1 n d 2 

A:=^ai(s,x)— + - ^ <J l3 (s,x) 



i=l 



dxidxj 
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The matrix a(s,x) = (cr^x), 1 < i, j < n), is defined in the following way 

cr(s, x) = B T (s, x)B(s, x), B(s, x) := (bki(s, x), 1 < k,i < n). 

According to the Lemma 1 and definition of function r] e (x) for proof of 
invariance of set K it suffices to prove the next inequality 

Arj e {x)>0, xe{K 3e \K e ), s>0. 
It is not difficult to check the following properties of function uj e {x — z) 



d d d 2 d 2 

(u € (x - z)) = — (cj £ (x - z)); - (u £ (x - z)) = - — — {io e (x - z)); 



r\ \ w fc\"^ 7 7 O V c V 7 7" O O V fc V 7 7 O O 

OZi OXi OZiOZj OXiOXj 

d d 2 

LO e (x - Z)\\ x _ z \ =e = Q^(u e (x - Z ))\\ x _ z \ =e = Q Xi Q ("e{x ~ z))\\x-z\=e = 0- 



We define set 



K e (x) :={z: \z - x\ < e n (K 3e \ K e )} 



Further, applying the properties of function uj e {x) and Taylor-series 
expansion of functions a^(s, x), <r^(s, x) in point z we get 



Me (X) = CLi(s ,X) — + - ^ <Tij(8 , x) ^ ^ J / X (z)uJ £ (x - z)dz 



/ x i \ k 



1 ^ d a,i(s,0(x,z)) , 3 , 



k,j 
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4^ fe£ l (it " zt)(im " zm>+ 

fc,m 

1 ^ d 3 (T ij (s,9 1 (x,z)) . u u \<9 2 w e (:r-z) 

+ - > « [Xk - z k ){x m - z m ){xi - zi) — — — dz 

6 f—' i dz k dz m ozi / dz t dzj 

k,m,l J J 

f f ® t t \ t w / .daAs,z) 

= ~ j ^^-^{^{S^^eix- Z))-U e (x- Z)—^ + 

K e {x) 1 



,xr d fr \ ( ,da,i(s,0(x,z))\ , . 

+ Z^~Q^. y- Xk ~ z k) u t\ x ~ z ) J + U A X - z ) 



a>i(s,z) 



- ^(x k - z k )u e (x - z) dz _ dzk + 

k,j 



d(Tij{s,z) du e (x - z) ^ d_ ( _ dojji^z) du e (x - z) \ 

+ f)y. \ ^ Xk Ay, f)y. + 



dzi dzj dzi \ " ' dzk dzj 



d<Tij(s, z) duj e {x - z) ^ _ d 2 aij(s,z) du e {x - z) 
dzi dzj ^ dzidzk dzj 
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+ ((x k -z k )(x -z ) ^i2^A diQe(x-z) \ 

+ 2^dz i [} k k)[ m m) dz k dz m dz 3 ) + 

u\ST(rr v ,d 2 a XJ (s,z)dLj e (x-z) | 



2/ .\ III "'llLf r\ r\ r\ 

OZiOZ m OZj 
m J 



L-X^(<r ~ , d 2 <Ji 3 (s,z) du € (x- Z) 
' 2 ^ " Zk) dz k dzi ' 3z 3 



-S^(r,-zuMr -z \ ^i2^A d ^(x-z) 

9 ^ Xk Z k)K X m Z rn ) p. p. p. o 



^— ' OZiOZkOZ m OZn 

k,m J 

d 3 a ij (s,9 1 {x,z)) .d 2 u e (x-z] 



k,m,l 



k,m 

+ E dxkdXmdxi i*« - **)(»m - *»)(«! - *») dzidzj }dz. (7) 



Here 6(x, z) = z + 6(x - z), 6i(x, z) = z + 6 x (x - z)\ < 9 U 9 < 1. 

Applying Ostrogradskii - Gauss theorem we transform some integrals in right 
part of (7) to integrals on the surface A e (x). Further, we obtain the estimate 
of smallness for some surface integrals and the volume integrals. 

We set 



,/• \ 2e 2 (^ - Zj) du e {x - z) 

f e {i,x,z) = 7^ : now = f e {i,x,z)uj e {x - z). 

{e z — \x — z\ z y ozi 

So 



/ XI ^7 ( a i( s > z )ve(x ~ z))dz = J ^2ai(s,z)vi(z)u e (x - z)d(3 z . (8) 

K e (x) 1 ' A £ (x) i 

We make use of Cuachy- unyakovskii inequality in (9) and later on. 

dai(s,z) 



I / (j2(x k ~ z k )u e (x - z)' 

K e {x) i * ^ k 



dz k 
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A e (x) 



< 



A e (x) 



\ 



EE"*) 



k \ i 



daj{s, z) 



uj e {x — z)d(3 z < eci, 



C\ < oo. 



(9) 



E E(^ - **>* - z > dz \ dZk dz \ * 

t \ i k 



K e (x) 



< 



^2(xk - z k y 



K e (x) 



E 



C 2 < OO. 



d 2 Oi(0(x,z)) du e (x-z) 

-{x k - z k ){xj - Zj) — dz\ < 



(10) 



dz k dzo 



dzi 



x ^y^J*(i,x,z)u e (x - z)dz < ec 3 , c 3 < oo. 



(11) 



We exploit the condition (2) for estimate in (12) 



A e (x) 
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^1 J ^^b ki (s,z)vi(z)^b kj (s,z)^d(3 z \ < 



A e (x) k i 



< - 



a; ( x) * p V i 

xcj e (x - z)d(3 z = o(l). (12) 
We exploit the relation cr^-(s,z) = aji(s,z) for estimate in (13) 



2 y , dzidzk dzj 

K e (x) 



1,3 fe J 



+ ly^, _^ ^ 2 cr z j(g,z) - z) 

2 m m dzidz m dzj 

m J 



l r 5 2 (j^(s,z) du e {x - z)\ 

k ' J 



lr^, . d 2 (T i j(s i Z) \ 



2 ^ ctefcctei 



1 /" s ^( d 2 <j ij {s,z) l d 2 <jjj{s,z) 1 d 2 cr zj (s,z) \ _ 

+ 2 J M a^a*,- " 2 a^a^- 2 a^a^ 7^ zj + 

K e (x) hj 

+1 / 3 a 3 ^-(5,z) i ^(s,z) l av^g^u ^ ^ l 

2 7 ^ V dzjdzidzk 2 dzjdzidzk 2 dzjdzkdzi J ~\ 

K e (x) 
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A e (x) V j V k 



£ 



dzidz k ) 



u e (x — z)dz-\- 



+ 



K e (x) 



— < C4€, C4 < OO. 



(13) 



d^aij{s, z) duJejx - z) 

dZkOZm dZj 



\ 



dz\ = 



/ 



A e (a:) 



j fc,m 



dzkdz 



xu e (x — z 



W,\<\ J £ (x* - **) 2 /£y?(j.*.*) /£^, 2 Wx 

A e (a) fc V J V * 



\ . ^ V dz k dz m 

\ i,J,k,m 



d(3 z < €€5, C5 < 00. 



(14) 



K e (:r) ' 



< 



/ <9 3 <7ij(s, 



y . 

\ \dzidz k dz 

\ J,k,m 



UJ e {x — z)dz < €Cq, C5 < OO. (15) 
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We exploit the following relation 

d 2 uj e {x — z) 



dzidzj 



= f e (i,X,z)f e (j,X,z) ~ 



2e 2 8 
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8e 2 (xi - Zi){xj - Zj) 



(e 2 — \x — z\ 2 ) (e 2 — \x — z\ 2 Y 



u e (x - z), 



here 5ij be Kronecker's symbol, 
for estimate in (16). 



1 1 f d 3 a ij {s,9 1 (x,z)) d 2 u e (x-z) 

12 I J E dzkdzkdzi -*»)(**-*») a a <**l< 



-12 / (X^-- z -) 2 ) \ ( 

X e (a;) ^ m ' ' i,j,k,m,l 



d 3 a ij (s,9i(x : z)) 
dz k dz m dzi 



X 



2e 2 V^ 8e 2 E(^"^) 2 



(e 2 — |x — z| 2 ) (e 2 — |x — z| 2 )' 



Cq < OO. 

Concider remaining summand in (7). 



(16) 



We 



d 



K e (x) 



1 (xk — Zk) — — 

<9z z V dz k 



daij(s J z) du e (x — z) 



dbpAz) 



d 



A, (.r) k V J " i 
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A, (.,•} k p j i 

= :h+h. 

The first summand ii is estimated analogy to (11) with help condition 2 of 
theorem. Thus we have I\ = o(l). 

We observe that by construction the points of boundary dA e (x) of set A e (x) 
have the following properties: 

z G dA e (x) =^> |x — z\ = e =4> w e (x — z) = 

Now we make use of local property of surfase dK e for more precise 
representation of summand I2 ■ 

The variables Zi, i = l,n in A e (x) have form Zi = yi, % < n — l,z n = 
Put y = (2/1, ._. . , 2/„), here y n = <£>(yi, . . . , y n _i). 

The domain D e {x) which corresponds to A e (x) has the following form 

D e (x) = {(yi, ■ ■ ■ , yn-i) ■ \x-y\< e}. (17) 

The boundary dD e (x) is set of points for which in (17) the next equality is 
fulfilled. Let y' = (yi, . . . , y n _i). Thus if y' G dD e (x), then w e (x — y) = 0. 
Set w e (x — y) = 0, under y' ^ D e (x). 

Thus the function u; e (iE — y) is finite function in space i? n_1 with 
support D e (x). The following formula of integration by parts is true for such 
functions [8, p. 106]. 

/q f d 

f—u e (x-y)dy' = - J u e (x - y)—fdy' z = l,n-l. 

D e (x) D e (x) 

(18) 

Applying (18) to integration in I2, we get 
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D e (x) j k P > 1 



— '■ 121 + ^22- 



To estimate of summand I21 with help Cauchy - Bunyakovskii's inequality 
we will use the condition 1 and the supposition that surface belongs to class 
C 2 

Later on it is convenient to omit the argument of functions. 



/21I < 



EE 



D e (x) 



J V,i 



E ^7 ( h PO y i^- h Pi I ( X k ~ Vk) 



d_ 



u e (x - y)dy' < 



< 



J to*-*?) '??(?{( 



D e (x) 



j p,i \ k 



dyj 1 dyk 



dvi db i d 2 b i ] \ 2 

+ ( & pj^^ 1 ) 2 + (fryj^ Q r ) 2 \ u e (x - y)dy < ec 2 i. 
PJ dyj dy k PJ dy 3 dy k \) 

Here C21 is bounded constant. 

Combining (8)-(16) and latter one gives the following representation 



Ar] e (x) = - J lY^ai(s,y)i/i(y)- 

D e (x) ^ 1 
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~\ £ b iv{^V>i{v)^f-^.Vk \ - y)dy' + "(I)- (19) 

j,k,p,i J 

Let G(s,y) denote the second summ in braces of right part of latter equality. 
It takes place the following relation for new variables 



dy n ~ fri ^ dyi 



(20) 



It is not hard to calculate the following equalities for partial derivatives in 
summands from G(s,y) 



f d 



db p i d 



dyk dy. 



Vk 



dy k u P*i 

0, 

d l 
db 



Xj 5 



dyk ^ Xk ' 



n—1 
d l 2 

m=l 



if j = k < n 

if j ^ k, j < n, k < n 
if j < n, k = n 

ii k < n, j = n 
ii j = k = n. 



Now we will show that the function G(s, y) coincides with the following 
function from condition 3 of theorem completely 



2^ — is t {z)b kj {s,z), s > 0; 



dzj 

i,j,k 

in case when differentiation is fulfilled in coordinates y. 

Applying (20), we get the following equailities for differentiation in (21) 



(21) 



db 



pi 



if j < n 



dzn 



n — 1 



n-l 
/ = 1 



It is clear that latter one defines the summands in (21) which is identical to 
the summands in the G(s,y). 
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Thus it follows from representation (19) that under conditions theorem there 
exists such e* > that the inequality Arj e (x) > is fulfilled for all e < e*. 
Theorem is prooved. 
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